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Abstract. The coupled Maxwell-Lorentz system describes feed-back action of electro- 
magnetic fields in classical electrodynamics. When applied to point-charge sources (viewed 
as limiting cases of charged fluids) the resulting nonlinear weakly hyperbolic system lies 
beyond the scope of classical distribution theory. Using regularized derivatives in the 
framework of Colombeau algebras of generalized functions we analyze a two-dimensional 
analogue of the Maxwell-Lorentz system. After establishing existence and uniqueness of 
solutions in this setting we derive some results on distributional limits of solutions with 
5-like initial values. 

Keywords: Regularized derivatives, Colombeau algebras, self-interaction, weakly hyper- 
bolic, Maxwell-Lorentz 

AMS subject classification: Primary 35Dxx, 46F10 , secondary 35Lxx, 78A25 
1. Introduction 

In recent years there has been considerable development in the application of nonlinear 
theories of generalized functions to problems in physical field theory, most successfully 
in general relativity (cf. 0, fL4|, |I5| , |17| ) but also starting in electrodynamics (see 



||10||). Due to the nonlinearities inherent in the coupling of fields and equations of motion 
methods from classical distribution theory soon turn out to be insufficient for treating 
such problems. On the other hand, the theory of algebras of generalized functions in 
the sense of J.F. Colombeau offers the possibility of applying large classes of nonlinear 
operations to distributional objects. 

We shortly recall the principal setup of electrodynamic field theory. Configuration space is 
modeled by 4-dimensional Minkowski space. The basic physical quantities are represented 
by the charge- and current density Jq, J, the electric field E, and the magnetic field B. 
Maxwell's field equations are 

(HI) div£ = (H2) d t B + curl E = 
(II) div£ = J (12) d t E — curl 5 = — J 

with initial conditions 

E\ t=0 = E B\ t=0 = B . 
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These equations imply the following well-known compatibility relations 



(CO) divB = 0, div£ = J \ 
(CI) d t J + dw.J = 0. 



t=Q 



On the other hand (CO), (CI) together with (H2), (12), and the initial conditions are 
sufficient to derive equations (HI), (II) (take the divergence of equations (H2), (12) and 
integrate with respect to time). So, in order to solve the complete Maxwell system it is 
sufficient to study the hyperbolic Cauchy problem consisting of initial data and equations 
(H2), (12) while E Q , B , J, and J are assumed to obey the relations (CO), (CI). 

Let us now investigate the nature and role of the charge density Jo and the current J. In 
an idealized model we will interpret these as quantities of a charged fluid (or sometimes 
clouds of charged dust, cf. |12|]) with charge density proportional to mass density (for 



convenience we set this ratio equal to 1). In this model we write 

(Jo, J) = p(uq,u) , (1.1) 

where p represents charge density distribution in space-time and (uo,u) is the 4-velocity 
field of the fluid, normalized with respect to the Minkowski metric by 

u 2 -\u\ 2 = l, (1.2) 

which is plausible as the 4-velocity is always assumed to be time-like. 

If we study the feed-back action of the electromagnetic field on the charged fluid the 
equations of motion for u are determined by the Lorentz force according to (cf.[|P2j|) 

3 

d t u + V ; Uj djU = E + - U x B . 



Observe that the above formulae involve products of the basic physical quantities (resp. 
of their coefficients). These quantities can be singular if for example jump discontinuities 
of certain conductor-dielectrica configurations are to be modeled or if point charges are 
studied. 

For instance, if one is interested in the equations of motion of charged particles resulting 
from the electromagnetic field they generate then p - and consequently E and B - involve 
strong singularities corresponding to concentration of charge on the particles' world line. 

Therefore in considering the realistic problem of self-interaction we have to deal with non- 
linear differential equations involving singular functions and products and compositions 
thereof. 

Even worse, the equations of motion cannot be written down consistently in the frame- 
work of classical distribution theory. Therefore, often physical motivations are given for 
renormalization and adding corrective terms. Unfortunately these "tricks" give rise to 
strange behaviour of solutions to the resulting equations (e.g. the Lorentz-Dirac equa- 
tion). Critical discussions of the emerging inconsistencies and the resulting unphysical 
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solutions ( "run-away-solution" , "preaccelaration" , . . . ) can be found in the books [[LJJ 
and ill. 



The aim of this paper is to analyze the mathematical structure of the coupled Maxwell- 
Lorentz system by means of a low dimensional "toy model" which displays some of the 
main features of the general problem. The plan of exposition is as follows: In section 2 we 
develop a (l+l)-dimensional model of self-interacting charges which will provide the basis 
of our further investigations. As the model equations are nonlinear, of evolution type and 
involve singular initial data we use the formalism of regularized derivatives in Colombeau 
algebras of generalized functions (cf. J7j). Section 3 presents the main existence and 
uniqueness result in this framework. Finally, in section 4 we consider singular initial 
data (point charges). We illustrate the influence of the regularization procedure on the 
existence of distributional limits of the unique solutions established above. 

2. A (l+l)-dimensional model of self-interacting charges 

We try to mimic the purely mathematical structure of Maxwell equations with only one 
space dimension, i.e. in a two-dimensional Minkowski space M with Lorentz metric in 
coordinates (t,x) equal to dt 2 — dx 2 . The field quantities (electric and magnetic) are 
scalar functions on M and will be denoted by E and B. Further "physical" quantities 
are represented by the scalar charge density p and the relativistic (normalized) velocity 
field (a/1 + m 2 , u) of a "fluid" consisting of a collection of (strictly time-like world lines of) 
charged particles. Therefore p(a/1 + u 2 , u) represents the relativistic current of the fluid. 
Interpreting divergence as simple ^-derivative and removing curl terms in the original 
Maxwell equations this produces the electromagnetic fields according to 

d x B = d t B = 



d t E = —pu d x E = p\J\ + u 2 . 

Therefore we may set B equal to some constant B . The second line defines a gradient 
condition on E which is solvable if and only if the following (relativistic) charge conser- 
vation holds 

d t (pVl + v?) + d x (pu) = . (2.1) 

In order to recover an evolution type system we simply combine the equations for E into 
one by adding them. 

Finally the feed-back action of the field on the charged particles of the fluid is modeled 
by mimicking the structure of the Lorentz force law: 

d t u+ , =d x u = E+ - B. (2.2) 

The left hand side evaluated along a world line z with z(s) = u(z(s)) gives z(s), i.e. an 
actual acceleration. 

Again, the above formulae involve products (and even more complicated nonlinear oper- 
ations) of the basic physical quantities. By the transformation p t— * o = pyl + u 2 we can 
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rewrite the system of equations in evolution form 











u 


H 









(i+« 2 ) 3 / 2 VTT^ 

Recall that E, u, and a are real valued functions of the variables (t, x) and B Q is a real 
constant. We study the Cauchy problem for system (|2.3|) with initial data 









u 


H 




-J 




V ^0 



when t = 0. (2.4) 



System ( |2.3|) is quasi-linear hyperbolic but not strictly hyperbolic. Also, its coefficients 
are not uniformly bounded. Furthermore, the main application we have in mind is to 
model self-interaction of charged particles (approximated by fluid clusters shrinking to 
isolated points) which qualitatively corresponds to gq being proportional to a delta dis- 
tribution. In order to address such questions we need a mathematical theory which is 
capable of describing singular objects and allows for an unrestricted application of both 
differentiation and of a sufficiently large class of nonlinear operations. Differential alge- 
bras of generalized functions containing the space of distributions as a subspace and the 
space of smooth functions as a faithful subalgebra and providing the above-mentioned 
tools were introduced by J.F. Colombeau in |J (see also [f|, []lT|). We will therefore work 
in this framework. 

Using the relation d x (l + w 2 ) 1 / 2 = u(l + m 2 ) -1 ' 2 ^^ and denoting a(y) = y(l + y 2 )^ 1 ^ 2 we 
can rewrite the above system in the more compact form 

d t E + d x E = a(l-a{u)) (2.5) 
d t u + d x (l + m 2 )5 = E + B a{u) (2.6) 
d t a + d x (aa(u)) = 0. (2.7) 



This form of the equations will be our starting point for the method of regularized deriva- 
tives which is the subject of the following section. 



3. Regularized derivatives — existence and uniqueness of solutions 
in Colombeau algebras 

We will work with a variant of Colombeau algebras denoted by (? Sjg (f2) (Q sm open subset 
of R n ) which is a simplified version with representatives globally bounded and smooth (cf. 
0): denote by r>Loo(f2) the algebra of restrictions of smooth functions on R n to Q which 
have all derivatives in L°°(R' 1 ); then within the algebra £ S)g (f2) = I\° o (^) ( ' ' oo ' ) we define 

£u,s,M = {(ue)e e £ S ,M | Va G No 3p > : \\d a u £ \\ hoo(n) = 0(e- p ) (e -> 0)} 
K,M = e ^(H) | Va G N^Vg > : \\d a u £ \\ Loom = 0{e q ) (e -> 0)} ; 



4 



then Q Sig (Q) is the quotient algebra of £M,s,g(^) niodulo Af Stg (Q). In our application we 
will have Q = (— T, T) xS for T > arbitrary. 

The notion of regularized derivative was introduced in 0-0. It extends weak solution 
concepts in algebras of generalized functions. The idea is to not only represent singular 
functions by (classes of) regularizations but also to consider differentiation as a limiting 
case of more smooth operations, namely convolution with derivatives of delta regulariza- 
tions. This means that a partial derivative djU of the generalized function U = cl[(u £ ) £ ] 
will be replaced by cl[(<9j^> e * u E ) e ) where (ip e )e is a delta net (in view of application to 
our model only first order derivatives are considered here). Thus if U is associated to a 
distribution T the regularized derivative will be associated to dfT. 

Regularized derivatives provide a powerful solution concept for a wide variety of partial 
differential equations. For example, the Cauchy problem for quasilinear evolution type 
systems 

d t u(t,x)= Aa(u(t,x))&*u(t,x)+B(u(t,x)) 

\a\<m 

u(0, x) = Uq(x) 

is uniquely solvable in this framework provided that A a and B are bounded in all deriva- 
tives (0, Thm. 5.1). 

The choice of which of the occurring derivatives are to be regularized may depend on 
special information about the model, the structure and form of the equations, and the 
type of effects one wants to observe. In our case we will systematically replace all x- 
derivatives in the model equations of the form (|2.5[ - |2.7|) but leave the time derivatives 



unchanged in order to reflect the evolution type of the equations. Note that this is 
different from regularizing x-derivatives in ( |2.3|) . 

The regularized x-derivative is defined by choosing a function ^eD with integral equal to 
1 on R - we call this a mollifier henceforth - and an increasing function h: (0, oo) — > (0, oo) 
with h(e) —* as e — ► 0. Setting (p v (x) = tp(x/v)/v we have d x (cp„(x)) = (f' v (x) = 
<p' (x I v) I v 1 (by (p' v we will always understand {tp v )') and define the regularized derivative 
of a generalized function U = cl[(u e ) £ ] by the formula 

(d x ) h U = d[{<p' Ke) *u e ) e ]. (3.1) 

The function h can be used to control the speed of convergence towards the usual derivative 
compared to the regularization parameter e. That it is well-defined and also the basic 
properties of the regularized derivative follow from the estimates 

\\d a d x ip h{£) *u £ \\ Loo < ||^|| Ll ||5 a S a: 'u £ || Loo (3.2) 

\\d x <Ph(e) * U e \\ hOB < ^(^^IIvIIliIKIIl- ■ ( 3 - 3 ) 

Now we can state the main existence and uniqueness result for the Cauchy problem with 
regularized x-derivatives. We will only have to impose the following growth condition on 
h: 
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\/p G N Bk G N s.t. exp(h(e)- p ) = 0(e- k ) (e -> 0) (*) 
For example, h(e) = CTn(|ln(e)|) -1 satisfies (*). 

Theorem 1. Let T > and assume that h satisfies (*). Then given Eq, uq, ctq in t? Sjg (R) 
the system 

d t E + x ) h E = a(l-a(u)) (3.4) 
d t u + (d x ) h (l + u 2 )^ = E + B a(u) (3.5) 
d t a + {d x ) h {aa{u)) = 0, (3.6) 

with initial conditions 

E | t=0 = Eq, u | t=0 = uq, a | t=0 = oq (3.7) 
has a unique solution (E,u,a) in Q S;g ([—T,T] x M) 3 . 

Proof. We will proceed in two steps: First, the candidate for a Colombeau solution will be 
built up by nets of smooth solutions at fixed e. Second, existence and uniqueness in ^ Sig will 
be proved by verifying additional estimates with respect to the regularization-parameter 
e. To begin with, we consider global solvability of the corresponding integro-differential 
equations where we will suppress all references to the index e. 

Step 1: integrating with respect to t and inserting initial conditions we get an equivalent 
system of integral equations 

E \ */ E(T,.)*tf/ \ 

«o p)- / (l + u(T,.) 2 ) 1/2 *<p' \(x)dr + 
cr J { \a(u(T,.))a(r,.) * <p' J 

* I <t(t,x)(1 - a(u(r,x))) \ 
+ / E(t,x) + B cl(u(t,x)) \dr. (3.8) 
I \ o J 

Writing V = (^1,^2,^3) instead of (E,u,o~) and denoting 

||V|| T = max sup \vj(t,x)\ 

7=1,2,3 \ t \<T,x<=R 

and Vq = (Eq,u ,o~o), which is globally bounded by assumption, we consider the right 
hand side as an operator R on the set 

B T (V ) = {V G C([-T,T] x R) 3 I || V - V Q \\ T < 1} . (3.9) 

Using ||a|| LOO < 1, \/l + s 2 < 1 + s for s > 0, and Young's inequality for the convolutions, 
one derives from (13.81) 




\(RV) 3 (t,x) - (V ) 3 (x)\ < Ty\\ u (l + \\V\\ T )+T(2\\V\\ T + \B \) (j = 1,2,3) 
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which shows that RV G Bt{Vq) for T small enough. To see that R is actually a contraction 
on Bt(Vq) for small T we can estimate the terms occurring in \\RV — RW\\ T for V, 
W G B t (Vq) after applying Young's inequality and/or taking sup for all integrands. First, 
in the integral involving convolutions the non-trivial factors are handled as follows: by 
the mean value theorem applied to the function f(y) = a/1 + y 2 and ||/'|| Loo = ||a|| L °° < 1 
we have sup \y/T+vl- a/TT^H < \\V- W\\ T ; since also a' is globally bounded by 1 
and by definition of B t (Vq) we further estimate 

\a{v 2 )v 3 - a(w 2 )w 3 \ < \a(v 2 )(v 3 -w 3 )\ + \{a{v 2 ) - a(w 2 ))w 3 \ 
< \v 3 - w 3 \ + \v 2 - w 3 \\w 3 \ < \\V - W\\ T (1 + \\V \\ hoo ) ; 

For dominating the sup of the integrands in the second integral we simply observe 

\v 3 (l - a(v 2 )) - w 3 (l - a(w 2 ))\ < \v 3 - w 3 \ + \a(v 2 )v 3 - a(w 2 )w 3 \ 
\vi + B a(v 2 ) —w 1 — B a(w 2 )\ < \v x - w x \ + \B \\a(v 2 ) - a(w 2 )\ 

to get terms we already estimated above. Altogether for T small, depending on \B \, 
||Vo|| L oo, and ||y'|| L i the map R: Bt{Vq) — > Bt{Vq) is a contraction. Therefore a solution 
of the fixed point equation exists for small T. 

Assuming that V and W are two fixed points for T > arbitrary one can estimate their 
difference similar to the above except that only sup^jj is to be considered inside integrals 
with respect to r. A standard Gronwall argument then yields global uniqueness of the 
solution. An a priori estimate for a solution V in [— T, T] x R is derived by using the fixed 
point formula (|3.8|) , |(1 + s 2 ) 1 / 2 — 1| < \s\, taking sup x . inside the integrals, and again 
applying Gronwall's lemma: 

\\V\\ T < (II V \\ LOO + T(|H| L1 + |5 |))e^'^ +1 ) . 

Therefore existence of a global solution follows. Clearly this solution is bounded on 
[-T, T] x R by definition of B T (V ). 

Finally to get unique solutions in T'l oo ([— T, T] x M.) from initial data in T> T oo ( IK. ) WG 
proceed inductively. Set B^iVo) = Bt(Vq) with ||V^|| T0 = ||^|| T ; then for k G N we define 

B*(v ) = {v e c k ([-T,T] xr) 3 | ||y-y || T>fc < 1} 

where 

\\y\\ Ttk = \\y\\ T , k -i + nwdtvwr^ + p.vw^) . 

Then estimates can be done essentially by the same arguments as above because the 
nonlinear combinations of unknowns remain of the same type (e.g. all derivatives of the 
function a are again bounded and therefore application of Young's inequality, mean value 
theorem, and Gronwall's lemma are still possible). The weight factor T k for derivatives of 
order k in the above norms serves to control inner derivatives and assures again e.g. the 
contraction property for small T. This concludes the purely classical base for the next 
step. 
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Step 2: Existence and uniqueness of solutions in Q s , s ([— T, T] x R) 3 . 
In equation ( |3.8|) all involved functions except (p now carry an index e and (p is to be 
replaced by <fh(e)- We assume ((E 0e , u 0e , cr 0e )) e G £M jSig (IR) 3 . In order to establish existence 
of a solution we will show that the net (V s ) = ((E s , u £ , cr £ )) £ belongs to ^M,s, g ([ — T, T] xK) 3 . 

From the integral equations for (V e ) and using Gronwall-type arguments as above we 
conclude 

max sup \vi >£ (t, x)\ < 

1<«<3 \t\<T,x£R 

< C(\\V 0}E \\ LOO ,T, ||^|| Ll )| J Bo, £ |^(^~ 1 exp(3T/ i ( £ )- 1 ||^ll L i), 

which (by (*)) gives the ^M-estimates of order 0. The ^M-estimates for the x- derivatives 
of higher order follow inductively by similar arguments if in each step we write the x- 
derivative of any convolution product of with any function / as y}( m+1 ) * / and 
use (*). The ^M-estimates for the t-derivatives follow directly by induction from the 
differential equations. 

Concerning uniqueness, suppose that W = cl[(W £ ) £ ] is another solution in Q S;S ([—T,T] x 
R) 3 . This means that there exists some N = cl[(iV e )J G Af([-T,T\ x R) 3 with 

/ v i, e - w h£ \ / (f 1,5 - w 1;£ ) ^ \ 

d t v 2>£ - w 2>£ = - [(1 - uyi - (1 - wl £ f>] * ip' h{£) + 

\V 3 ,6 ~ Wz,e ) \ [a(v 2 , e )v 3t e ~ a(w 2 ,e)w 3 ,e) ) 

v 3>£ (l - o(«2, e )) - w 3 , e (l - a(v 3>e )) \ I N h£ \ 
vi,s ~ w h£ + B 0j e(a(v 2t£ ) - a{w 2>£ )) + N 2 , £ 
/ \N 3 J 

From this by arguments similar to those used above we inductively derive AA-estimates 
for each derivative of (V £ - W £ ) £ . Thus V = W in G s , g ([-T, T] x R) 3 . 

□ 



Remark 1. Having obtained a unique solution to the field equations let us reinvestigate 
the behaviour of point particles subjected to the velocity field (vT+ u 2 , u), i.e. determine 
its integral curves in R 2 . Classically, this amounts to considering the system 

z (s) = ^l + u(z(s)) 2 

zi(s) = u{z(s)) (3.10) 
2(0) = (t ,x ) 

where z — (zo, z±) : / — > R 2 (I some interval containing 0). By the first line of ( |3.10| ) we 
may reparametrize via r = zq(s), thereby obtaining the following equation for (r,w(r)) = 
z(s): 

w(r) = a(u(r,w(r))) (3-H) 
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Of particular interest to us is u arising as a solution of system ( p.4|) - (|3.6| ) with strongly 
singular (e.g. 5-like) initial data. In this case, (|3.11| ) is an example of a nonlinear ODE 



involving generalized functions. Note that its right-hand side contains a composition of 
generalized functions (defined by componentwise insertion of representatives). A theory of 
such equations in the framework of Colombeau generalized functions has been developed 
by R. Hermann and M. Oberguggenberger in ||, |9|]. Since all derivatives of a are globally 
bounded and u G Q Sjg it follows that cl[(a o u £ ) £ ] is in Q Stg and therefore of L°°-type (see 
||, Def. 2.3). Thus ||, Thm. 3.3 ensures the existence of a solution in Q which is 
actually in Q s g (which follows inductively from (|3.11|) ). In particular, this allows to derive 
(generalized) trajectories of charged particles without having to resort to renormalization 
procedures altering the actual equations of motion. This may indicate a way of avoiding 
the inconsistencies mentioned at the end of section 1. 

4. Point charges as initial data 

As we mentioned in the introduction we will now study the situation where the initial data 
(Eq, uq, ctq) are qualitatively of the form (0, 0, 5). This can be modeled by choosing cr to 



be the class of a strict delta net (p e ) £ in ^ S)g (M) in the sense of [[Tl|], Def.7.1. According 
to this definition (p e ) £> o is a family of test functions in D(IR) such that 

supp p £ — > {0} as e — > 
fp £ (x)dx = l for all e > (4.1) 
J\p £ {x)\ dx is bounded independently of e . 

By the above theorem we obtain a unique solution in Q s ,g([—T, T] x M) 3 of system ( |3.4l|3.6| ) 
with (|3.7p for appropriate scaling function h. However, the dependence of the definition 
of the regularized derivative on the mollifier (p is essential. This was observed already 
in || and illustrated by numerical tests. Also in the context of scalar conservation laws 
it was shown how the choice of "left" or "right" mollifiers allows to recover the classical 
entropic solutions. 

We are not able to give a complete description of conditions for the existence of distribu- 
tional limits of solutions of the more complicated system (|3.4| - |3.6| ). Nevertheless we will 
present some results analyzing the consequences of employing one-sided mollifiers. 

Proposition 1. Let T > 0, h as in Thm.l, and (E ,uo,ao) = (0, 0, cl[(p e ) £ ]) with 
(p £ ) £ a strict delta net. If suppy? C [— oo,0] (resp. C [0,oo]) then the unique solution 
(E,u,a) G <? s ,g([ — T, T] x E) to ( j3^4"l - [3.6|) with initial conditions ( |3.7|) is supported in 
[— T, T] x (-oo,0] (resp. [— T, T] x [0,oo)). 

Proof. It will suffice to consider the case of a left-sided mollifier (since a similar argument 
applies to the case of a right-sided mollifier). Fix some xq > 0. By our assumption on the 
support of <f we have 

* o 

a £ (t, x) = - —-r J J (p'(y)a{u £ (r, x - h(£)y))a £ (r, x - h(£)y)dydr 
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for x > xq and e sufficiently small (depending on xq). Also, since a is globally bounded 

t 

C f 

sup \a £ (t } x)\ < ——\\(p'\\ L1 / sup \a £ (t,x)\dr, 

X>XO "H^V J X>XQ 



so sup x>;! , \(t e (t, x) \ = by Gronwall's inequality which (since x was arbitrary) establishes 
the claim for a. Then since 

t 

C f 

sup \E £ (t, x)\ < sup \a e (t,x)(l — a(u e (t,x)))\ + —, J I^'Hl 1 / SU P \E e (t, x)\dr, 

x>xq x>xq llyE J J x>xo 



it follows that sup x>xo \E e (t,x)\ = 0. Finally, 

t o 

f 



d x u £ {t, x) = - — r / / (p' (y)a(u £ (T, x - h(e)y))d x u £ (r, x - h(e)y)dydr + 



o -1 

t 



+ B a'(u £ (T,x))d x u £ (r,x)dr 



o 



From this it follows that d x u £ (t,x) = for x > 0. But then since (1 + u £ ) ' * V'hU) 
(a(u £ )d x u £ ) * <Ph(e) — in this region we obtain 

t 

d t u e (t,x) = B a(u £ (t,x)) = B / a'(u e (r, x))d t u e (r, x)dr . 



o 



Thus also d t u e (t, x) = and therefore u £ (t, x) = for x > 0. □ 

Applying this result we now demonstrate how the choice of a left mollifier can destroy the 
intuitive expectation of field propagation along the characteristics of the original system. 



Corollary 1. The (unique) solution of ( |Q| - |3.6D , ( p.7| ) according to Prop.l cannot be 
associated to a distribution in D'(R 2 ) 3 . 

Proof. Assume that (E, u, a) £ D'(IR) 3 is the distributional limit of (E £ , u £ , a £ ). Then from 
d t E £ + (p'h{£) * E £ = a £ (l — a{u £ )) it follows that a £ a{u £ ) — > a — E, where F = d t E + d x E. 
Also, the third line of the regularized system implies d t a + d x (a — F) = 0, so 

d t (a - d x E) + d x (a - d x E) =0 (4.2) 

Using the coordinate transformation T(t, x) = (x,t — x) it follows that a — d x E is of the 
form v(t-x) (i.e. T*(l<g>t>)) for some v G V'(R), hence is contained in C°°(R, V'(R)). By 
Prop.l the supports of a and E are contained in x < and from the initial conditions we 
get (a — d x E)(0) = 8(x). Thus a — d x E = 5(t — x), which contradicts supp(er — d x E) C 
{x < 0}. □ 

It remains an open question whether there are choices of the mollifier (p which yield 
distributional limits of the solution for singular initial data as in Prop.l. 
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Remark 2. Let us interpret equations ( |2.5[ - |2.7|) in the form P(E, u, cr) = with the 
operator P : C°°(M 2 ) 3 -> C°°(M 2 ) 3 given by 

/ d t E + d x E + a(a(u) - 1) \ 
P(£, u, cr) = <9 4 w + a(u)d x u -E- B a(u) . (4.3) 
\ d t a - d x (a(u)a) ) 

We can decompose P into a sum of a linear differential operator L and a (purely) nonlinear 
part Q with 

' + d x E - a ' 
L(E,u,a) =| j (1.1) 

d t cr 

a(u)a 

Q(E,u,a) = | a(u)(d x u-B ) J . (1.5) 
-^(a^cr) 



We observe that in a linearization of P (in the sense of [ 13| , p. 20) at a point (toj^o) 
with w(t ,^o) = & H terms corresponding to Q vanish and we are left exactly with 
L. Therefore, if we consider some smooth approximation of the initial data (with q the 
constant point charge) 

(E ,u ,a ) = (0,0, q5) (4.6) 

and a linearization of the system near t = 0, the unique distributional solution (E, u, a) £ 
C 1 ([0,T\,V'(R)) 3 of the Cauchy problem 



(4.7) 



with initial condition ( |4.6|) may give a hint on how to construct approximate smooth 
solutions of the non-linear system. The solution to the linearized system is given by 

H(x) - H(x-t) \ 
(t-x)(H(x) -H(x-t)) , (4.8) 
5(x) J 

where H denotes the Heaviside function. However, using a standard regularization of this 
solution (by means of a model delta net) as a candidate for such an approximate solution 
runs into difficulties: explicit calculations show that the blow up of the term aa(u) at 
x = is unavoidable unless regularizations with 5H « are considered. 









: 


(t,x)=q- 
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